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Abstract 

We report a first-principles study on the geometric, electronic and magnetic properties of V- 
doped anatase TiC>2. The DFT+U (Hubbard coefficient) approach predicts semiconductor band 
structures for Tii_ x V x 02 (x=6.25 and 12.5 %), in good agreement with the poor conductivity of 
samples, while the standard calculation within generalized gradient approximation fails. Theo- 
retical results show that V atoms tend to stay close and result in strong ferromagnetism through 
superexchange interactions. Oxygen vacancy induced magnetic polaron could produce long-range 
ferromagnetic interaction between largely separated magnetic impurities. The experimentally ob- 
served ferromagnetism in V-doped anatase TiC>2 at room temperature may originate from a combi- 
nation of short-range superexchange coupling and long-range bound magnetic polaron percolation. 
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The discovery of ferromagnetism (FM) in Co-containing Ti02 semiconductors has at- 
tracted significant attention because of its potential application for developing functional 
dev.ces that .nan.pulate both 8 pi„ and cha.ge.fl However, the origin of FM observed in 
Co-doped anatase Ti02 at and above room temperature is still under debate. I2| Many ex- 

nnn 

periments argued the intrinsic nature of FM, y, \4, M while the presence of the Co clusters in 

nil 

the samples could not be completely excluded. |6|, |7[ To study the observed FM, theoretical 



efforts have been made through analyzing the electronic and magnetic structures for vari- 
ous Co configurations in the host matrix. Yang et al. have pointed out that ferromagnetic 
ordering occurs only when Co atoms locate close to each other. 8J One plane- wave pseudopo- 
tential calculation supports that the short-range exchange between adjacent Co atoms could 



contribute to FM 
the FM coupling. [10 




n_addition, oxygen vacancy is shown to play a key role in enhancing 
12] Recently, V-doped anatase Ti02 has been discovered to have 
a surprisingly high Curie temperature {Tc > 400 K). 13| In contrast to Co-doped Ti02 
that has been widely studied theoretically, few investigations on V-doped Ti02 have been 
reported in the literature. One interesting work by Wang et al. using LDA and LDA+U 
methods investigates the electronic structure of Tii_ x V r 02 (x=6.25%) and argues that the 
neighboring V-V cations prefers FM at a high doping concentration (x=25%) by employing 
a 2x1x1 supercell. [l^] However, the relative stability among magnetic phases is extremely 
sensitive to atomistic structure, its associated defects' distribution and dopant concentra- 
tion, further investigation is desirable to gain better understanding of high Tc FM in this 
system. 

In this paper, we use ab initio band-structure and total energy methods, implemented in 
the Vienna Ab initio Simulation Package (VASP), 15|, [l6| to study the geometric, electronic 
and magnetic properties of V-doped anatase Ti02 ( Tii_ x V x 02), in particular the effects 
of doping concentration (x=6.25% and 12.5%) as well as oxygen vacancy. The projector 
augmented wave (PAW) method is chosen to represent the electron-ion interaction. 15J] Ex- 
change correlation interactions are described by the Perdew-Burke-Ernzerhof generalized 

n 

gradient approximation (GGA).[17[ The Brillouin-zone (BZ) integration is performed on a 



well converged Monkhost-Pack k-points grid 



The plane wave kinetic energy cutoff is set 



to be 400 eV. Atomic positions and lattice parameters are optimized at GGA level until the 
atomic forces are smaller than 0.01 eV/A. To account for the strongly correlated interactions 
of the V 3d electrons, a moderate on-site Coulomb repulsion (U=3.0 eV) has been applied 
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only to V 3d orbitals since further correction for the host material has little impact on the 
magnetic ordering. 

We start with a Ti 15 V0 32 supercell to model the low doping concentration case (6.25%), 
where a 48-atom 2x2x1 supercell of pure anatase Ti02 (Ti 16 032) is demonstrated in Fig.l 
and one V atom substitutes Ti at Ti\ site. The lattice parameters of optimized Ti 15 V032 
supercell are a=7.64 A and c=9.65 A, which are slightly smaller than those of perfect TiC>2 
(our theoretical lattice parameters of Ti^C^ are a=7.65 A and c=9.68 A at GGA level). 
Clearly, this type of substitution only leads to a small geometrical distortion in the vicinity 
of V impurity (less than 0.1 A). To discuss the solubility of V in anatase Ti02 host, it 
is useful to define the substitution energy of V impurity as E s =E(Ti 15 V0 32 )+E(Ti)-E(V)- 
E(Ti 16 32 ), where E(Ti 16 32 ),E(Ti) and E(V) is the total energy of a 2x2x1 supercell of 
pure anatase Ti02, the bulk hep Ti and bec V, respectively. E s is predicted to be 2.34 eV, 
which is even less than one third of the substitution energy of Co dopant (8.51 eV). 8| This 
result agrees nicely with experimental observation that V impurity is well dissolved into 
Ti0 2 .|l3| 

The band structure, total density of states (DOS) and V 3d partial DOS (PDOS) of 
Ti 15 V0 3 2 are presented in FigJ21 The calculation at GGA level suggests that the DOS near 
the Fermi surface {Ep) mainly originates from V t2 9 orbitals, which is shown in Figj2](c) 
for clarity. The V impurity band has strong anisotropic character in the a-b plane (i.e. 
A^R and T— >M), however, it is less dispersive along the c axis of the BZ [i.e. M— sA and 
Z— >r). It is important to note that Ep is largely crossed by spin- up states, while spin-down 
states only cut Ep slightly. This indicates that the system is nearly half-metallic at GGA 
level. As there exists strongly correlated interaction in V 3d shells, we choose DFT+U as a 
scheme beyond-GGA level. The calculated PDOS by DFT+U are plotted in Figs. 2(e) and 
2(f). We find that the metal- insulator-transition occurs when the parameter U is applied 
to the V atom. This strong correlation lifts the degeneracy of t2 9 states, so as to open the 
gap. An occupied band composed of mainly d xy component is split off from the bottom 
of the conduction band. These features suggest that Hubbard U correction increases the 
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ization of the V d electrons. Results presented above are fundamentally consistent with 



Jjj] except for their lifted Ti PDOS due to additional corrections for Ti 3d states. 
It can be expected that the direct impurity interaction can not yield experimentally ob- 
served strong FM because of the large V-V distance. The original supercell is extended 
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along b and c axis to form 2x4x1 and 2x2x2 supercells (Ti3oV2064) to examine magnetic 
ordering. From Table I, the calculated energy difference between the AFM and FM states 
(AE) is negative, which indicates that the ground state is AFM for both cases as expected. 
Moreover, the on-site Coulomb repulsion at V atom strengthens AFM with a small increase 
of magnetic moment. The AFM ground state is further confirmed by our calculation using 
hybrid B3LYP functional as implemented in CRYSTAL03 code. 20J To understand the ex- 
perimental observed high Tq in V-doped TiC>2, one has to search other origins to lead FM. 
In the rest of this paper, first we explore the effects of closely distributed V impurities esp. 
at high doping level, then oxygen vacancy on electronic and magnetic properties. 

According to Janisch et a/.'s suggestion, we have constructed three Ti^V^C^ struc- 
tures to model high impurity concentration (12.5%) of V-Ti02, where a second V atom 
occupies Ti 2 , Ti 3 or Ti 4 site in the TiisVC^ supercell, named the chain, grid and nearest- 
neighbor (NN) configurations (the optimized V-V distance is 3.82, 5.39 and 3.06 A), respec- 
tively. From Table I, AE is positive for all three cases within GGA approach and increases 
with V-V distance decreasing. We find that the energy of the NN configuration is the lowest, 
strongly favoring V-dopants clustering. 

The DOS shapes have almost the same characteristic shapes among three configurations, 
hence, only the one for the nearest-neighbor one is presented in Figj3l Similar to the Co-Ti02 
case,[9j when x increases from 6.25% to 12.5%, the total DOS, V 3d PDOS and magnetic 
moment of V cation change slightly. The DOS and PDOS by DFT+U are presented in 
Figs. 3(b) and 3(d), which reveal that the system becomes semiconductor when the parameter 
U is applied. It is consistent with the poor conductivity observed in V-Ti02 samples. The on- 
site U correction shifts d xy state down to the middle gap, while at GGA level it is a mixture 
of three t2 9 orbitals (d xy , d yz and d xz ) locating at Ep. The larger separation between the 
occupied and empty V 3d states, due to strong correlation, results in the decrease of the AE 
value and even the change of its sign for the chain and grid configurations in Table I . The 
remarkable observation in our work is the high fidelity of AE for the NN configuration. As 
it is extremely hard to choose the accurate value for U, this result is particularly valuable 
to highlight the important contribution of V-O-V NN configuration to yield FM state. 

The reason for the weakened FM couplings within DFT+U could be sought in terms 
of the p-d hopping mechanism [l^]. The V impurities interact with orbitals of the same 
symmetry and form a set of bonding-antibonding states. Within GGA approximation, FM 
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state has an energy gain through the partial occupancy of V £25 states. The Hubbard U 
at V sites results in a fully occupied d xy band. In a FM arrangement, both bonding and 
antibonding levels are filled, and there is no energy gain for this coupling. In the AFM 
arrangement, the bonding states are filled for both spin up and down channels, while the 
antibonding ones are empty. It implies that AFM is favored between adjacent V cations 
when correlation effect is taken into account. However, there exists a strong competition 
due to superexchange. It is well-known that superexchange may have a FM contribution 
for filled shells. Recently, Janisch et al. argue that short-range superexchange between 90° 
metal-anion-metal bond plays a significant role in transition metal doped anatase Ti02- 9j 
We believe the FM originates from the superexchange interaction with the V-O-V bond 
angle about 104.9° in the NN configuration. 



As indicated by previous work 
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121 ]. the presence of oxygen vacancy (Vo) may influ- 



ence the distribution of dopants. To study the impact of Vo, three cases are considered: 
(a) Ti 16 03i with one oxygen atom removing from the 2x2x1 supercell of TiC^; (b) The 
basal oxygen (Oi) of V-contained octahedron is removed from TiisVC^ and subsequently 
the oxygen deficient cell is doubly enlarged along b axis, named as the basal Oi 2x4x1 con- 
figuration; (c) Vo occupies the vertex site of V-incorporated octahedron (by removing O2) 
and the original supercell is doubly extended along c axis, denoted as the vertex O2 2x2x2 
configuration. After removing an oxygen atom from the cation-contained octahedron, the 
original structure experiences a considerable distortion. Cations next to the Vo site are 
repelled away. The calculated total energies indicate that Vo prefers to stay close to V than 
near Ti. 

Figj4]shows the DOS results for the Ti 16 03i and basal Oi 2x4x1 configurations. When 
an oxygen atom is removed, a nonmagnetic defect state (mainly from Ti cations near Vo) 
appears at the bottom of the conduction band. If a V atom occupies the cation site of 
the oxygen-deficient octahedron, V 3d states overlap in energy with the defect state which 
yields 0.11 e transferred from the defect state to the empty V 3d spin-up state enhancing 
spin polarization of V ion. It fits quite well with the physical picture proposed by Coey 
et al. that the hybridization between the defect state and magnetic dopants' states will 
promote FM.j^l] However, the Vo mediated magnetic coupling is sensitive to distance. 22] 
The polarons induced by Vo tend to parallelize their spins, corresponding to a ferromagnetic 
arrangement, only when their distance is less than a critical value, which is the case in the 
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basal Oi 2x4x1 configuration. When their distance increases larger than this critical value, 
the strict spin alignment between polarons disappeared. This leads to a loss of long-range 
coherence, as in the vertex O2 2x2x2 case. Further efforts are needed to determine the 
exact critical distance required for polaron percolation in this system. DFT+U calculations 
qualitatively do not change the result given by GGA. 

In summary, we have studied the structural, electronic and magnetic properties of V- 
doped anatase Ti02 by both GGA and DFT+U calculations. The correlation effect is not 
well represented by GGA, while DFT+U calculations give a better physical picture. The 
results suggest that direct exchange between well separated V ions cannot account for the 
observed FM in this system. The superexchange interaction is responsible for the stable 
FM between nearest-neighbor V cations. The long-range interactions between oxygen de- 
fects induced magnetic polarons are sensitive to their distance. Provided that a macroscopic 
ferromagnetic clusters appear by polarons overlapping, the system favors FM. Thus, a com- 
bination of the long-range polaron percolation and the short-range superexchange through 
nearest-neighbor V atoms could contribute to the high Tq FM in V-doped anatase Ti02- 
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TABLE I: Calculated energy difference between the AFM and FM states (AE=E(AFM)-E(FM)), 
and the ground state magnetic moment (M) per V atom for various configurations at GGA level 
and at DFT+U level. The data inside the parentheses are obtained by DFT+U approach. 

Unit cells AE(meV/atom) M(fM B ) 



T130V2O64: 

2x2x2 -1.0 (-4.2) tO. 71 (±0.83) 

2x4x1 -8.0 (-12.7) ±0.69 (±0.80) 
Tii4V2032: 

chain 60.2 (-5.0) 0.84 (±0.81) 

grid 21.8 (-56.0) 0.77 (±0.83) 

nearest-neighbor 75.2 (14.5) 0.89 (0.91) 
T130V2O62: 

basal Oi 2 x 4 x 1 18.5 (18.5) 1.14 (1.37) 

vertex 2 2 x 2 x 2 -3.8 (-4.5) ±0.97 (±1.17) 
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FIG. 1: (Color online) The schematic 48-atom supercell of anatase TiO"2 (2x2x1). The gray and 
red balls stand for Ti and O atoms, respectively. 
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FIG. 2: (Color online) Band and DOS structure of V-doped anatase TiO"2 modeled by a 2 x 2 x 1 
supercell (TiisVC^). (a),(b) the majority and minority spin GGA band structure, (c) majority 
spin around E F . Total DOS and V 3d PDOS (d) GGA results, (e) DFT+U results, (f) DFT+U 
results of the projected DOS of V 3d orbitals. The Fermi energy is shifted to zero. r=(0,0,0), 
M=(l/2,l/2,0), A=(l/2,l/2,l/2), R=(l/2,0,l/2), and Z=(0,0,l/2). 
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FIG. 3: (Color online) DOS structure of the nearest-neighbor configuration (Tii4V2 032) with a 
ferromagnetic arrangement of the V moments. Total DOS and V 3d PDOS (a) GGA result, (b) 
DFT+U result. Projected DOS of V 3d orbitals (c) GGA result, (d) DFT+U result. The vertical 
line denotes the Fermi level. 
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FIG. 4: (Color online) DOS structure of oxygen deficient TiigOsi and the basal Oi 2 x 4 x 1 
configuration (Ti3oV2062)- (a) Total DOS and 3d PDOS of Ti cations near the oxygen defect for 
Tii 6 3 i. Total DOS and V 3d PDOS for the basal Oi 2 x 4 x 1 obtained within (b) GGA approach, 
(c) DFT+U method. The vertical line denotes the Fermi level. 
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